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Abstract 

A lattice Boltzmann model is developed by coupling the density (D2Q9) and the temperature distribution functions 
with 9-speed to simulate the convection heat transfer utilizing Al 2 0 3 -water nanofluids in a square cavity. This 
model is validated by comparing numerical simulation and experimental results over a wide range of Rayleigh 
numbers. Numerical results show a satisfactory agreement between them. The effects of Rayleigh number and 
nanoparticle volume fraction on natural convection heat transfer of nanofluid are investigated in this study. 
Numerical results indicate that the flow and heat transfer characteristics of Al 2 0 3 -water nanofluid in the square 
cavity are more sensitive to viscosity than to thermal conductivity. 



List of symbols 

c Reference lattice velocity 
c s Lattice sound velocity 
c p Specific heat capacity (J/kg K) 
e a Lattice velocity vector 
f a Density distribution function 

f£1 Local equilibrium density distribution function 

F a External force in direction of lattice velocity 
g Gravitational acceleration (m/s ) 
G Effective external force 
k Thermal conductivity coefficient (Wm/K) 
L Dimensionless characteristic length of the square 
cavity 
Ma Mach number 
Pr Prandtl number 
r Position vector 
Ra Rayleigh number 
t Time (s) 

T a Temperature distribution function 

T eq Local equilibrium temperature distribution 

function 

T Dimensionless temperature 

To Dimensionless average temperature (T 0 = (T H + T C )I2) 
T H Dimensionless hot temperature 
T c Dimensionless cold temperature 
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u Dimensionless macrovelocity 

w c Dimensionless characteristic velocity of natural 
convection 
w a Weight coefficient 
x, y Dimensionless coordinates 

Greek symbols 

P Thermal expansion coefficient (K~ ) 
p Density (kg/m 3 ) 

v Kinematic viscosity coefficient (m 2 /s) 
X Thermal diffusion coefficient (m 2 /s) 
H Kinematic viscosity (Ns/m 2 ) 
(p Nanoparticle volume fraction 
S x Lattice step 
S t Time step t 

Tf Dimensionless collision-relaxation time for the flow 
field 

t T Dimensionless collision-relaxation time for the tem- 
perature field 

AT Dimensionless temperature difference (AT = T H - 
T c ) 

Erro^ Maximal relative error of velocities between 
two adjacent time layers 

Error 2 Maximal relative error of temperatures between 
two adjacent time layers 

Subscripts 

a Lattice velocity direction 
avg Average 
C Cold 
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f Fluid 
H Hot 
nf Nanofluid 
p Particle 

Introduction 

The most common fluids such as water, oil, and ethylene- 
glycol mixture have a primary limitation in enhancing the 
performance of conventional heat transfer due to low ther- 
mal conductivities. Nanofluids, using nanoscale particles 
dispersed in a base fluid, are proposed to overcome this 
drawback. Nanotechnology has been widely studied in 
recent years. Wang and Fan [1] reviewed the nanofluid 
research in the last 10 years. Choi and Eastman [2] are the 
first author to have proposed the term nanofluids to refer 
to the fluids with suspended nanoparticles. Yang and Liu 
[3] prepared a kind of functionalized nanofluid with a 
method of surface functionalization of silica nanoparticles, 
and this nanofluid with functionalized nanoparticles have 
merits including long-term stability and good dispersing. 
Pinilla et al. [4] used a plasma-gas-condensation-type clus- 
ter deposition apparatus to produce nanometer size- 
selected Cu clusters in a size range of 1-5 nm. With this 
method, it is possible to produce nanoparticles with a 
strict control on size by controlling the experimental con- 
ditions. Using the covalent interaction between the fatty 
acid-binding domains of BSA molecule with stearic acid- 
capped nanoparticles, Bora and Deb [5] proposed a novel 
bioconjugate of stearic acid-capped maghemite nanoparti- 
cle with BSA molecule, which will give a huge boost to the 
development of non-toxic iron oxide nanoparticles using 
BSA as a biocompatible passivating agent. Wang et al. [6] 
showed the method of synthesizing stimuli-responsive 
magnetic nanoparticles and analyzed the influence of glu- 
tathione concentration on its cleavage efficiency. Huang 
and Wang [7] produced e-Fe 3 N-magnetic fluid by chemi- 
cal reaction of iron carbonyl and ammonia gas. Guo et al. 
[8] investigated the thermal transport properties of the 
homogeneous and stable magnetic nanofluids containing 
y-Fe203 nanoparticles. 

Many experiments and common numerical simulation 
methods have been carried out to investigate the nano- 
fluids. Teng et al. [9] examined the influence of weight 
fraction, temperature, and particle size on the thermal 
conductivity ratio of alumina-water nanofluids. Nada et al. 
[10] investigated the heat transfer enhancement in a hori- 
zontal annuli of nanofluid containing various volume frac- 
tions of Cu, Ag, A1 2 0 3 , and Ti0 2 nanoparticles. Jou and 
Tzeng [11] studied the natural convection heat transfer 
enhancements of nanofluid containing various volume 
fractions, Grashof numbers, and aspect ratios in a two- 
dimensional enclosure. Heris et al. [12] investigated 
experimentally the laminar flow-forced convection heat 
transfer of Al 2 0 3 -water nanofluid inside a circular tube 



with a constant wall temperature. Ghasemi and Aminossa- 
dati [13] showed the numerical study on natural convec- 
tion heat transfer of CuO-water nanofluid in an inclined 
enclosure. Hwang et al. [14] theoretically investigated the 
natural convection thermal characteristics of Al 2 0 3 -water 
nanofluid in a rectangular cavity heated from below. 
Tiwari and Das [15] numerically investigated the behavior 
of Cu-water nanofluids inside a two-sided lid-driven differ- 
entially heated square cavity and analyzed the convective 
recirculation and flow processes induced by the nanofluid. 
Putra et al. [16] investigated the natural convection heat 
transfer characteristics of CuO-water nanofluids inside a 
horizontal cylinder heated and cooled from both of ends, 
respectively. Bianco et al. [17] showed the developing lami- 
nar forced convection flow of a water- A1 2 0 3 nanofluid in a 
circular tube with a constant and uniform heat flux at the 
wall. Polidori et al. [18] investigated the flow and heat 
transfer of Al 2 0 3 -water nanofluids under a laminar-free 
convection condition. It has been found that two factors, 
thermal conductivity and viscosity, play a key role on the 
heat transfer behavior. Oztop and Nada [19] investigated 
the heat transfer and fluid flow characteristic of different 
types of nanoparticles in a partially heated enclosure. Ho 
et al. [20] carried out an experimental study to show the 
natural convection heat transfer of Al 2 0 3 -water nanofluids 
in square enclosures of different sizes. 

The lattice Boltzmann method applied to investigate the 
nanofluid flow and heat transfer characteristic has been 
studied in recent years. Hao and Cheng [21] simulated 
water invasion in an initially gas-filled gas diffusion layer 
using lattice Boltzmann method to investigate the effect of 
wettability on water transport dynamics in gas diffusion 
layer. Xuan and Yao [22] developed a lattice Boltzmann 
model to simulate flow and energy transport processes 
inside the nanofluids. Xuan et al. [23] also proposed 
another lattice Boltzmann model by considering the exter- 
nal and internal forces acting on the suspended nanoparti- 
cles as well as mechanical and thermal interactions among 
the nanoparticles and fluid particles. Arcidiacono and 
Mantzaras [24] developed a lattice Boltzmann model for 
simulating finite-rate catalytic surface chemistry. Barrios 
et al. [25] analyzed natural convective flows in two dimen- 
sions using the lattice Boltzmann equation method. Peng 
et al. [26] proposed a simplified thermal energy distribu- 
tion model whose numerical results have a good agree- 
ment with the original thermal energy distribution model. 
He et al. [27] proposed a novel lattice Boltzmann thermal 
model to study thermo-hydrodynamics in incompressible 
limit by introducing an internal energy density distribution 
function to simulate the temperature field. 

In this study, a lattice Boltzmann model is developed by 
coupling the density (D2Q9) and the temperature distribu- 
tion functions with 9-speed to simulate the convection 
heat transfer utilizing nanofluids in a square cavity. 
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Lattice Boltzmann method 

In this study, the Al 2 03-water nanofluid of single phase 
is considered. The macroscopic density and velocity 
fields are still simulated using the density distribution 
function. 



where c$ = is the lattice sound velocity, and w a i, 

p h a is the weight coefficient. 

The macroscopic temperature field is simulated using 
the temperature distribution function: 



/.(r + . A ,t + a,)-/.(r,0 = ~[/.(r.O-/."'('.0] + *A (D T a (r + e a S t ,t + S, ) - T a ( r,t) = -±-[ T „ ( r,t)-T? (r, t)] (6) 



K -») f eq 



(2) 



where if is the dimensionless collision-relaxation time 
for the flow field; e a is the lattice velocity vector; the 
subscript a represents the lattice velocity direction; f a (r, 
t) is the population of the nanofluid with velocity e a 

(along the direction a) at lattice r and time t; ( r, t ) 

is the local equilibrium distribution function; S t is the 
time step t; F a is the external force term in the direction 
of lattice velocity; G = -P{T nC -T 0 )g is the effective exter- 
nal force, where g is the gravity acceleration; f} is the 
thermal expansion coefficient; T is the temperature of 
nanofluid; and T 0 is the mean value of the high and low 
temperatures of the walls. 

For the two-dimensional 9-velocity LB model (D2Q9) 
considered herein, the discrete velocity set for each 
component CC is 



(0,0) 



,[(a-l)|],rin[(a-l)| 



4lc\ cos| (2a -l)— .sin (2a -1) — 



a =0 

a = 1, 2, 3, 4 
a =5,6,7,8 



(3) 



where c = 8* / 8 t is the reference lattice velocity, 8 X is 
the lattice step, and the order numbers a = 1, 4 and 
a = 5, 8, respectively, represent the rectangular direc- 
tions and the diagonal directions of a lattice. 

The density equilibrium distribution function is cho- 
sen as follows: 



1 + - 



u | (e a -uY 



2c; 



2c; 



(4) 



where t t is the dimensionless collision-relaxation time 
for the temperature field. 

The temperature equilibrium distribution function is 
chosen as follows: 



T a eq = w,T 



1 + 3 i_-^ + 4. 5 ^ X "J -1.5-" 



2c" 



2c 1 



(7) 



The macroscopic temperature, density, and velocity 
are, respectively, calculated as follows: 



a=0 



5> 

/ Ja e a 



(8) 



(9) 



(10) 



a=0 



The corresponding kinematic viscosity and thermal 
diffusion coefficients are, respectively, defined as follows: 



v = — c 
3 



Z= 3 C 



M - " \S t 



(ID 



(12) 



For natural convection, the important dimensionless 
parameters are Prandtl number Pr and Rayleigh number 
Ra defined by 



Pr 



(13) 



4 
9 
1 
9 
1 

36 



a = 0 

a = 1,...,4 

a = 5,... ,8 



(5) 



Ra 



gPATL Pr 



(14) 



where AT is the temperature difference between the 
high temperature wall and the low temperature wall, 
and L is the characteristic length of the square cavity. 
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Another dimensionless parameter Mach number Ma 
is defined by 



Ma ■ 



u c 



(15) 



where w c = J g/3ATL is the characteristic velocity of 

natural convection. For natural convection, the Boussi- 
nesq approximation is applied; to ensure that the code 
works in near incompressible regime, the characteristic 
velocity must be small compared with the fluid speed of 
sound. In this study, the characteristic velocity is 
selected as 0.1 times of speed of the sound. 

The dimensionless collision-relaxation times Tf and x T 
are, respectively, given as follows: 



• 0.5 + 



• 0.5 + 



MaLJ3Pr 



c^dt^Ra 

3v 
Prc 2 dt 



(16) 
(17) 



Lattice Boltzmann model for nanofluid 

The fluid in the enclosure is Al 2 0 3 -water nanofluid. 
Thermo-physical properties of water and A1 2 0 3 are given in 
Table 1. The nanofluid is assumed incompressible and no 
slip occurs between the two media, and it is idealized that 
the Al 2 03-water nanofluid is a single phase fluid. Hence, the 
equations of physical parameters of nanofluid are as follows: 
Density equation: 



Pnf = (l-0)Pf +« 



(18) 



where p nf is the density of nanofluid, 0 is the volume 
fraction of A1 2 0 3 nanoparticles, pbf is the density of 
water, and p p is the density of A1 2 0 3 nanoparticles. 

Heat capacity equation: 



c pnf = (l-0)Sf 



PP 



(19) 



where C pn f is the heat capacity of nanofluid, C p f is the 
heat capacity of water, and C pp is the heat capacity of 
A1 2 0 3 nanoparticles. 

Dynamic viscosity equation [28]: 



Pnf 



Pf 



{l-cp) 



2.5 



(20) 



Table 1 Thermo-physical properties of water and Al 2 0 3 [29] 



Physical properties 


Fluid phase (water) 


Nanoparticles (Al 2 0 3 ) 


p (kg/m 3 ) 


997.1 


3970 


c p (J/kg K) 


4179 


765 


M (m 2 /s) 


0.001004 


/ 


k (Wm/K) 


0.613 


25 



where fi n f is the viscosity of nanofluid, and pi( is the 
viscosity of water. 
Thermal conductivity equation [28]: 



fe nf = kf 



(fe p + 2fe f )-2y(fe f -fe p ) 
(fcp + 2fc f ) + <p(fcf-fc p ) 



(21) 



where k n f is the thermal conductivity of nanofluid, and 
kf is the thermal conductivity of water. 
The Nusselt number can be expressed as 



Nu 



hH 



The heat transfer coefficient is computed from 
u — «w 



(22) 



(23) 



The thermal conductivity of the nanofluid is defined 



by 



Kt = - 



dT /dx 



(24) 



Substituting Equations (23) and (24) into Equation 
(22), the local Nusselt number along the left wall can be 
written as 



Nu = - 



dT 



II 



dx ) T H -T L 
The average Nusselt number is determined from 



(25) 



1 



(26) 



Results and discussion 

The square cavity used in the simulation is shown in 
Figure 1. In the simulation, all the units are all lattice 
units. The height and the width of the enclosure are 
all given by L. The left wall is heated and maintained 
at a constant temperature {T H ) higher than the tem- 
perature (7c) of the right cold wall. The boundary 
conditions of the top and bottom walls are all adia- 
batic. The initialization conditions of the four walls are 
given as follows: 



= 0 u = 
0 u = 



0, T = 1; x = 
0, dT I dy = 0; 



1 U : 



0, 
1 u 



T = 0 

-- 0, dT I dy = 0 



(27) 



In the simulation, a non-equilibrium extrapolation 
scheme is adopted to deal with the boundary, and the 
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adiabatic surface 



v 



— 



adiabatic surface 



X 



Figure 1 Schematic of the square cavity 



standards of the program convergence for flow field and 
temperature field are respectively given as follows: 



Erroq 



^JMm,- + M-Mm,Q]MMm.' + M-",(W)] 1 } 



Error, = 



^ T x {i,j,t + S t ) 2 



h (28) 



(29) 



where e is a small number, for example, for Ra = 8 x 
10 4 , ej = 10. 7 , and e 2 = 10_ 7 ; for 7?a = 8 x 10 s , ej = 10_ 8 , 
and e 2 = 10_ 8 . 

In the lattice Boltzmann method, the time step t = 1.0, 
the lattice step S = 1.0, the total computational time of 
the numerical simulation is 100 s, and the data of equili- 
brium state is chosen in the simulation. 

As shown in Table 2, the grid independence test is 
performed using successively sized grids, 192 x 192, 256 
x 256, and 300 x 300 at Ra = 8 x 10 s , <p = 0.00 (water). 
From Table 2, it can be seen that the numerical results 
with grids 256 x 256 and 300 x 300 are more close to 
those in the literature [20] than with grid 192 x 192, 
and there is little change in the result as the grid 
changes from 256 x 256 to 300 x 300. In order to 



Table 2 Comparison of the mean Nusselt number with 
different grids 

Physical 192 x 192 256 x 256 300 x 300 Literature 
properties [20] 



accelerate the numerical simulation, a grid size of 256 x 
256 is chosen as the suitable one which can guarantee a 
grid-independent solution. 

To estimate the validity of above proposed lattice 
Boltzmann model for incompressible fluid, the model 
is also applied to a nanofluid with nanoparticle volume 
fraction (p = 0.00 in a square cavity, and the research 
object and conditions of numerical simulation are set 
the same as those proposed in the literature [20]. Fig- 
ure 2 compares the numerical results with the experi- 
mental ones, and a satisfactory agreement is obtained, 
which indicates that it is feasible to apply the model to 
incompressible liquids with good accuracy. In Figure 2, 
there are a few differences because the nanofluid in 
the simulation is supposed as a single phase, while the 
real nanofluid is a two-phase fluid. Therefore, the 
small differences are accepted in the simulation, and 
the model is appropriate for the simulation of 
nanofluid. 

Figure 3 illustrates the velocity vectors and isotherms 
of the Al 2 0 3 -water nanofluid at different Rayleigh num- 
bers with a certain volume fraction of A1 2 0 3 nanoparti- 
cles (<p = 0.00). It is observed that there are two big 
vortices in the square cavity at Ra = 8 x 10 s ; as the Ray- 
leigh number increases, they are less likely to be 
observed compared with the condition at smaller Ray- 
leigh numbers. This may be because of the gradually 
increasing Rayleigh number (corresponding to the 
increase of the velocity), which causes the nanofluid to 
rotate mainly around the inside wall of the square cav- 
ity. In addition, it can be seen that the temperature iso- 
therms become more and more crooked as Ra increases, 
which illustrates that the heat transfer characteristics 
transform from conduction to convection. 

Figures 4 and 5 present the velocity vectors and iso- 
therms at Ra = 8 x 10 and Ra = 8 x 10 for various 
volume fractions of A1 2 0 3 nanoparticles, respectively. 



Nl4, 



8.367 



7.915 



7.704 



16 
15 
14 
13 
12 
11 
6010 

I 9 
8 



Hoet al.[20] 
Present work 




1000000 2E6 3E6 4E6 5E6 

Ra 

Figure 2 Comparison of the mean Nusselt number at different 
Rayleigh numbers. 
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Figure 5 Velocity vectors (on the left, ->0.002) and isotherms (on the right) for Al 2 0 3 -water nanofluid at Ra = 8 x 10 with different 
volume fractions, (a) <j> = 0.00, (b) tj> = 0.01, (c) <j> = 0.03, (d) <j> = 0.05. 



There are no obvious differences for velocity vectors and 
isotherms with different volume fractions of nanoparti- 
cles, which is because the volume fractions are so small, 
it is not significant in this case on comparing with Ray- 
leigh number, and the effect of those volume fractions is 
negligible. However, it can be seen that there is a little 
difference on local part of the isotherms, for example, as 
the volume fraction of A1 2 0 3 nanoparticles increases, 
the lowest isotherm in Figure 4 and the second lowest 
isotherm in Figure 5 become less and less crooked, 
which indicates that high values of (f> cause the fluid to 
become more viscous which causes the velocity to 
decrease accordingly resulting in a reduced convection. 
It is more sensitive to the viscosity than to the thermal 
conductivity for nanofluids heat transfer in a square cav- 
ity. This phenomenon can also be observed in Figure 6. 

Figure 6 illustrates the relation between the average 
Nusselt number and the volume fraction of nanoparticles 
at two different Rayleigh numbers. It is observed that the 
average Nusselt number decreases with the increase of 
the volume fraction of nanoparticles for Ra = 8 x 10 and 
Ra = 8 x 10 s . In addition, it can be seen that the average 
Nusselt number decreases less at a low Rayleigh number. 
For the case of Ra = 8 x 10 4 and Ra = 8 x 10 5 , it is indi- 
cated that the high values of <f) cause the fluid to become 
more viscous which causes reduced convection effect 
accordingly resulting in a decreasing average Nusselt 
number, and the flow and heat transfer characteristics of 
nanofluids are more sensitive to the viscosity than to the 
thermal conductivity at a high Ra. 



Conclusion 

A lattice Boltzmann model for single phase fluids is 
developed by coupling the density and temperature dis- 
tribution functions. A satisfactory agreement between 
the numerical results and experimental results is 
observed. 

In addition, the heat transfer and flow characteristics 
of Al 2 0 3 -water nanofluid in a square cavity are investi- 
gated using the lattice Boltzmann model. It is found 
that the heat transfer characteristics transform from 
conduction to convection as the Rayleigh number 
increases, the average Nusselt number is reduced with 
increasing volume fraction of nanoparticles, especially at 



-»-Ra= 


ixlO 4 


-»-Ra= 


ixlO 5 




0.00 0.01 0.02 0.03 0.04 0.05 



Figure 6 Average Nusselt numbers at different Rayleigh 
numbers 
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a high Rayleigh number. The flow and heat transfer 
characteristics of Al 2 0 3 -water nanofluid in a square cav- 
ity are demonstrated to be more sensitive to viscosity 
than to thermal conductivity. 
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